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Abstract 

We have performed a detailed investigation of total lifetimes for the H+ + and 
baryons in the framework of operator product expansion over the inverse mass of charmed 
quark, whereas, to estimate matrix elements of operators obtained in OPE, some approx- 
imations of nonrelativistic QCD are used. This approach allows one to take into account 
the corrections to the spectator decays of c-quarks, which reflect the fact, that these 
quarks are bound, as well as the contributions, connected to the effects of both the Pauli 
interference for the H+ + -baryon and the weak scattering for the H+-baryon. The real- 
ization of such program leads to the following estimates for the total lifetimes of doubly 
charmed baryons: r~++ = 0.43 ± 0.1 ps and r~+ = 0.11 ± 0.01 ps. 
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1 Introduction 



A study on weak decays of doubly charmed baryons is of great interest because of two 

reasons. The first one is connected to the investigation on the basic properties of weak 

interactions at the fundamental level, including the precise determination of CKM matrix 

parameters. The second reason is related to the possibility to explore QCD as it is 

provided by the systems containing the heavy quarks. In the limit of a large scale given 

by the heavy quark mass, some aspects in the dynamics of strong interactions become 

simpler and one gets a possibility to draw definite model-independent predictions. Of 

course, both these topics appear in the analysis of weak decays for the doubly charmed 

baryons, whose dynamics is determined by an interplay between the strong and weak 

interactions. That is why these baryons are the attractive and reasonable subjects for 

the theoretical and experimental consideration. 

i—i (*) 

The doubly charmed Sc C -baryon represents an absolutely new type of objects in 
comparison with the ordinary baryons containing light quarks only. The basic state of 
such baryon is analogous to a (Qg)-meson, which contains a single heavy antiquark Q 
and a light quark q. In the doubly heavy baryon the role of heavy antiquark is played 
by the (cc)-diquark, which is in antitriplet color-state. It has a small size in comparison 
with the scale of the light quark confinement. Nevertheless the spectrum of (ccg)-system 
states has to differ essentially from the heavy meson spectra, because the composed (cc)- 
diquark has a set of the excited states (for example, 25 and 2P) in contrast to the heavy 
quark. The energy of diquark excitation is twice less than the excitation energy of light 
quark bound with the diquark. So, the representation on the compact diquark can be 
straightforwardly connected to the level structure of doubly heavy baryon. 

Naive estimates for the lifetimes of doubly charmed baryons were done by the authors 
slightly early A simple consideration of quark diagrams shows, that in the decay of 
E+ + -baryons, the Pauli interference for the decay products of charmed quark and the 
valent quark in the initial state takes place in an analogous way to the D + -meson decay. 
In the decay of H+ , the exchange by the VF-boson between the valence quarks plays 
an important role like in the decay of D°. These speculations and the presence of two 
charmed quarks in the initial state result in the following estimates for the lifetimes: 

r(S++) « ±t(D+) ~ 0.53 ps, 

r(3+) « ~r(£>°) ~ 0.21 ps. 

In this work we discuss the systematic approach to the evaluation of total lifetimes 
for the doubly charmed baryons on the basis of both the optical theorem for the inclusive 
decay width and the operator product expansion (OPE) for the transition currents in 
accordance with the consequent nonrelativistic expansion of hadronic matrix elements 
derived in OPE. Using OPE at the first step, we exploit the fact, that, due to the 
presence of heavy quarks in the initial state, the energy release in the decay of both 
quarks is large enough in comparison with the binding energy in the state. Thus, we 
can use the expansion over the ratio of these scales. Technically, this step repeats an 
analogous procedure for the inclusive decays of heavy-light mesons as it was reviewed in 
0. Exploring the nonrelativistic expansion of hadronic matrix elements at the second 
step, we use the approximation of nonrelativistic QCD [||, Q, which allows one to reduce 
the evaluation of matrix elements for the full QCD operators, corresponding to the 
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interaction of heavy quarks inside the diquark, to the expansion in powers of where 
p c = m c v c ~ 1 GeV is a typical momentum of the heavy quark inside the baryon. The 
same procedure for the matrix elements, determined by the strong interaction of heavy 
quarks with the light quark, leads to the expansion in powers of ®S D . 

This way, taking into account the radiation of hard gluons in these decays, leads 
to the expansion in powers of a s , v c = and . It is worth to note, that this 

expansion would be well defined, provided the expansion parameters to be small. In 
the c — > sud transition, the ratio of typical momentum for the heavy quark inside the 
hadron to the value of energy, released in the decay, is not so small. We would like 
also to stress the important roles, played by both the Pauli interference and the weak 
scattering, suppressed as — ^ with respect to the leading spectator contribution, but the 
former ones are enhanced by a numerical factor, caused by the ratio of two-particle 
and three-particle phase spaces J5|. Numerical estimates show that the value of these 
contributions is considerably large, and it is of the order of 40 — 140%. These effects 
take place in the different baryons, H+ + and S+, and, thus, they enhance the difference 
of lifetimes for these baryons. The final result for the total lifetimes of doubly charmed 
baryons is the following: 

r =+ + = 0.43 ± 0.1 ps, 

1 — 'cc 

r = + = 0.11 ± 0.01 ps. 



2 Operator product expansion 

Now let us start the description of our approach for the calculation of total lifetimes 
for the doubly charmed baryons. The optical theorem, taking into account the integral 
quark-hadron duality, allows us to relate the total decay width of the heavy quark with 

the imaginary part of its forward scattering amplitude. This relationship, applied to the 

(*) 

Eic -baryon total decay width T„(*), can be written down as: 



-■(*) 
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where the Hcc state in ([l]) has the ordinary relativistic normalization, (S 
and the transition operator T is determined by the expression 

T = 5m [ d A x {TH eff (x)H eff (0)}, 



(*)|— (*)\ _ npv 
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where H e ff is the standard effective hamiltonian, describing the low energy interactions 
of initial quarks with the decays products, so that 
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where f is the number of flavors. 

Assuming that the energy release in the heavy quark decay is large, we can perform 
the operator product expansion for the transition operator T in (|I|). In this way we find 
a series of local operators with increasing dimension over the energy scale, wherein the 
contributions to T„(*) are suppressed by the increasing inverse powers of the heavy quark 
masses. This formalism has already been applied to calculate the total decay rates for the 
hadrons, containing a single heavy quark || (for the most early work, having used similar 
methods, see also @,@]). Here we would like to stress that the expansion, applied in this 
paper, is simultaneously in the powers of the inverse heavy quark mass and the relative 
velocity of heavy quarks inside the hadron. Thus, the latter points to the difference from 
the description of both the heavy-light mesons (the expansion in powers of ° ) and 
the heavy-heavy mesons [|| (the expansion in powers of relative velocity of heavy quarks 
inside the hadron, where one can apply the scaling rules of nonrelativistic QCD Q). 

In this work we will extend this approach to the treatment of baryons, containing 
two heavy quarks. The operator product expansion applied has the form: 

T = dirfcc + -Lc 2 (fi)cga^G^c + -^O(l). (4) 

The leading contribution in OPE is determined by the operator cc, corresponding to 
the spectator decays of c-quarks. The use of the equation of motion for the heavy quark 
fields allows one to eliminate some redundant operators, so that no operators of dimension 
four contribute. There is a single operator of dimension five, Qgq = Qgo'^uG^Q. As 
we will show below, significant contributions come from the operators of dimension six 
Q2Q2q = QTqqT Q, representing the effects of Pauli interference and weak scattering for 
S+ + and E+,, correspondingly. Furthermore, there are also other operators of dimension 
six Qqiq = Qa^jiD^G^Q and Qq 2 q = QD^G^ U T V Q. In what follows, we do not 
calculate the corresponding coefficient functions for the latter two operators, so that the 
expansion is certainly complete up to the second order of only. 

Further, the different contributions to OPE are given by the following 

T~++ = T 35c + % pj, 

1 — 'cc 1 

= %5c + %,WS, 

^cc 1 

where the first terms account for the operators of dimension three O3Q and five Oqq , the 
second terms correspond to the effects of Pauli interference and weak scattering. The 
explicit formulae for these contributions have the following form: 

T 35c = 2 * (r CjSpec cc - %(2 + KqJPx + K 2c P 2 ]0 Gc ), (5) 

where r 0c = and K 0c = + K 2c = 2(C+ — Ci). This expression has been 

derived in || (see also jT0|), and it is also discussed in p. The phase space factors P{ 
look like H 0: 

P l = {l-y)\ P 2 = (l-yf, 

where v = 

m c 

Pc.spec denotes the contribution to the total decay width of the free decay for one of 
the two c-quarks, which is explicitly expressed below. 
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For the effects of Pauli interference and weak scattering, we find the following for- 
mulae: 

([( (1 ~/" )2 - (1 ~/" )3 )fe7a(l " 75) Ci )(g j7 Q (l " 75)^) + 
,(l-z_) 2 (l-z_) 3 w w „, . 

C 2 - - 1 3 ' )(ci7a75Ci)(g 3 -7 a (l - 75)^)] (6) 
[(C+ + C_) 2 + ~(1 - ^)(5C 2 + C* - 6C_C+)] + 

K (1 ~/~ )2 - (1 ~/- )3 )(Q 7a (l - 75 )c i )(? i 7 a (l " 75)Q) + 



2 4 
,(l-z_) 2 (l-2_) 3 



)(c i7c( 75C i )(9 i 7 a (l " 75)%)]^(5C 2 + CI - 6C_C+)), 



Tws = 2 ^ P \(l-z + ncl + Cl + \(l-kh){Cl-Cl)) 

(Q7a(l " 75)c i )fe-7 a ( 1 - 75)^) + (7) 

fc§ (C <2 _ C 2 )( - 7q(1 _ 75 ) Ci ) fe . 7 a (1 _ 75 ) g .)] ) 

where p+ = p c + P q , P- = Pc - P q and z± = fc = a s (p)/a s (m c ). 

In the numerical estimates for the evolution of coefficients C+ and C_ , we have taken 
into account the threshold effects, connected to the 6-quark, as well as the threshold 
effects, related to the c-quark mass in the Pauli interference and weak scattering. 

In expression (|5|), the scale p. is approximately equal to m c . For the Pauli interference 
and weak scattering, this scale was chosen in the way to obtain an agreement between the 
experimental differences in the lifetimes of A c , and Hj3-baryons and the theoretical 
predictions, based on the effects, mentioned above. This problem is discussed below. 
Anyway, the choice of these scales allows some variations, and a complete answer to this 
question requires calculations in the next order of perturbative theory. 

The contribution of the leading operator cc corresponds to the imaginary part of 
the diagram in Fig. 1, as it stands in expression @. The coefficient of cc can be 
obtained in the usual way by matching of the Fig. 1 diagram, corresponding to the 
leading term in expression (||), with the operator cc. This coefficient is equivalent to the 
free quark decay rate, and it is known in the next-to-leading logarithmic approximation 
of QCD [12, 13, 14, 15, 16], including the strange quark mass effects in the final state 



[16]. To calculate the next-to- leading logarithmic effects, the Wilson coefficients in the 
effective weak lagrangian are required at the next-to-leading accuracy, and the single 
gluon exchange corrections to the diagram in Fig.l must be considered. In our numerical 
estimates we use the expression for T spec , including the next-to-leading order corrections, 
s-quark mass effects in the final state, but we neglect the Cabibbo-suppressed decay 
channels for the c-quark. The bulky explicit expression for the spectator c-quark decay 
is placed in the Appendix. 

Similarly, the contribution by Oqq is obtained, when an external gluon line is at- 
tached to the inner quark lines in Fig. 1 in all possible ways. The corresponding co- 
efficients are known in the leading logarithmic approximation. Finally, the dimension 
six operators and their coefficients arise due to those contributions, wherein one of the 
internal u or d quark line is "cut" in the diagram of Fig. 1. The resulting graphs are 
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depicted in Figs. 2 and 3. These contributions correspond to the effects of Pauli inter- 
ference and weak scattering. We have calculated the expressions for these effects with 
account for both the s-quark mass in the final state and the logarithmic renormalization 
of effective electroweak lagrangian at low energies. 

Since the simultaneous account for the mass effects and low-energy logarithmic renor- 
malization of such contributions has been performed in this work for the first time, we 
would like to discuss this question in some details. 

The straightforward calculation of diagrams of Figs. 2 and 3 with the account for 
the s-quark mass yields the following expressions: 



2G 2 F 2 „ (l-z-) 3 a0 , (l-z_T {l-z„)\ p°Lpi. 
^T P - [{ ^2— 9 + ( ~ 3~ ^ 



[(C+ + C^) 2 (c ila {l - 75 )q)(ct(1 " 75)<&) + (8) 
(5C 2 + - 6C+C_ + C 2 )(c* 7a (l - 75)c,-)(w(l " 75)%)], 

9G 2 m 2 

Tws = -^^pl(l-z + ) 2 [(Cl + C£)(c i7a (l- 75)^)^(1 -75)^) + 

(C 2 + - C 2 )(c na (l - 75)^)^7/3(1 - 75)%)]. (9) 
For p + and p- we use their threshold values: 

H -), P-=Pc{l -), 

taking into account that the logarithmic renormalization of effective low-energy la- 
grangian has the following form 

Leff,io 9 = ^{^ + c^ + i(i-*i)(c^.-c^)](cr /i )(^cO + 

~(C 2 - C 2 )kkcT,d)(d^c) + 1 -{C\ - Cl)kk{k^ - iWY^c)^ - (10) 
~[(C + + C„) 2 + i(l - k^)(5C 2 + + C 2 - 6C + C_)](cE> + h W c)(uT»u) - 
h^(5C 2 + + C 2 _- 6C + C_)(c i r M c fc + \<klvTtsc k )(u k T*Ui) - 

o o 

\[{C+ - C-f + \(\ - fc*)(5C£ + C 2 + 6C + C_)](cT> + |^c)(5T"5) - 

h^(5C 2 + + c 2 _ + ec + c_)(c i r /iCit + ^ 7At7 5c fe )(s fc r^) - 

o o 

^5(fcT - 1)(5C 2 + Ci)(cr m ^ a c+ ^c^ 75 t a c)j^}, 

where = 7/Lt (l - 75 ), k = (a s (p)/a s (m c )) and j° = u 7/i t a u + d 7/ ,i a (i + s 7/ ,t a s is the 
color current of light quarks (t a = X a /2 being the color generators). Having performed 
the manipulations, we have obtained formulae @, (Q). 

Here we would like to make a note, concerning the terms of effective lagrangian, 
containing the color current of light quarks. In the analysis below, we have omitted 
these terms, because they contribute in the lagrangian with the strength factor k~9 — 1, 
whose numerical value is equal to 0.054 (see below). 

To calculate the contribution of semileptonic modes to the total decay width of Hc C - 
baryons (we have taken into account the electron and muon decay modes only) we use 
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the following expressions [|1C|] (see also |l 

T sl = 4r c ({l-8p + 8p 3 



p 4 -12p z lnp} + 
E c {5 - 24p + 24p 2 - 8p 3 + 3p 4 - 12p 2 In p} + 
K c {-6 + 32p - 24p 2 - 2p 4 + 24p 2 In/)} + 
G c {-2 + 16p - 16p 3 + 2p 4 + 24p 2 lnp}), 
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where F c = \V CS \ 2 G 2 FJ ^, p 
by the expressions: 

G c = 



Ill- 
Ill 2 



The quantities E c = K c + G c , K c and G c are given 



;(*) 



2m^ 



(«)>, 



gG a pa' 
Am 2 



a/3 



:(*) 



(«)>, 



(12) 



where the spinor field c v in the effective heavy quark theory is defined by the form: 



c (x) 



-im c v-x 



1 + 



iD 



2m r 



Cy(x) 



(13) 



Thus, we can see, that the evaluation of total lifetimes for the doubly charmed baryons 
is reduced to the problem of estimation for the matrix elements of operators, appearing 
in the above expressions, which is the topic of next section. 



3 Evaluation of matrix elements. 

Let us calculate the matrix elements for the operators, obtained as the result of OPE 
for the transitions under consideration. In general, it is a complicated nonperturbative 
problem, but, as we will see below, in our particular calculation we can get some reliable 
estimates for the matrix elements of required operators. 

Using the equation of motion for the heavy quarks, the local operator cc can be 
expanded in the following series over the powers of ^-: 

/s?(*)|6cp(*)\ -1 ( E cc\c[^D) 2 ~ ^aG)]c\E c * c } } norm 1 

\— cc \ cc \^cc /norm — J- 2m 2 ^m 3 

Thus, this evaluation can be reduced to the calculation of matrix elements for the fol- 
lowing operators: 

c(iD) 2 c, (^)caGc, cj a (l - j5)cqj a (l - 75)?, cj a -f 5 cqj a (1 - j 5 )q. 

The first operator corresponds to the time dilation, connected to the motion of heavy 
quarks inside the hadron, the second is related to the spin interaction of heavy quarks 
with the chromomagnetic field of light quark and the other heavy quark. Further, the 
third and fourth operators are the four-quark operators, representing the effects of Pauli 
interference and weak scattering. 

In the system, containing the nonrelativistic heavy quark, the quark-antiquark pairs 
with the same flavor can be produced with a negligible rate, since the energy greater than 
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rriQ is required. In this situation, it is useful to integrate out the small components of the 
heavy-quark spinor field and to present the result in terms of the two component spinor 
^q. Following this approach, we find that all contributions from virtualities greater 
than fi, where m c > n > m c v c , can be explicitly taken into account in the perturbative 
theory. This method is general and analogous to the effective heavy quark theory. So, 



cc 



^ 9 ^ c -^t(D 9 m c+ ... (i5) 

cga^G^c = -2*+g*B* c - — *+(Dg£)* c + ... (16) 

In these expressions we have omitted the term ^>~^a{gE x D)^ c , corresponding to the 
spin-orbital interaction, because it vanishes in the ground states of doubly charmed 
baryons. By definition, the two-component spinor fy c has the same normalization as Q, 

cfix^t^t = J d 3 xQ + Q. (17) 

Then, with the required accuracy, ^ c can be expressed through the big components of 
spinor Q 

Q = e -imt f<P\ (18) 



X 



due to the following formula 



(this can be checked with the use of the equation of motion). Let us note that the 
covariant derivative should be taken in the adjoint representation, when it acts on the 
chromoelectric field, 

(DE) = {dT a - gf abc T b A c )E a . (20) 

Radiative corrections modify the coefficients of the chromomagnetic term (crB) and "Dar- 
win" term in (]l5|). However, in the situation at hand, these effects can be consistently 
neglected. 

Now let us consider the significance of different contributions to the expansions in 



(llq ) and (16). Evaluating the contributions of chromomagnetic and "Darwin" terms, we 
have to take into account the interaction of heavy quark with the light quark as well as 
the interaction with the other heavy quark. In the first case, the procedure of calculation 
is analogous to that for the heavy-light mesons. So, the "Darwin" term is suppressed by 
a factor of C ^ D in comparison with the chromomagnetic term, and, thus, we neglect its 
contribution. In the second case, the analysis is analogous to that for the heavy-heavy 
mesons, so that we can use the scaling rules of nonrelativistic QCD Q. In this approach, 
the contributions of different operators can be estimated, using the following relations 
in Coulomb gauge: 



i 

^ c ~ (m c v c )^, D ~ m c v c , gE ~ m 2 c vl, gB ~ m^vj, g ~ v? . 

From these scaling rules for the heavy-heavy interaction, we can deduce that the contri- 
bution of the "Darwin" term has the same order as that of chromomagnetic term. 
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Let us now start the calculation of matrix elements with the use of potential models 
for the bound states of hadrons. While estimating the matrix element value of the kinetic 
energy, we note, that the heavy quark kinetic energy consists of two parts: the kinetic 
energy of the heavy quark motion inside the diquark and the kinetic energy, related to the 
diquark motion inside the hadron. According to the phenomenology of meson potential 
models, in the range of average distances between the quarks: 0.1 fm < r < 1 fm, the 
average kinetic energy of quarks is constant and independent of both the quark flavors, 
constituting meson, and the quantum numbers, describing the excitations of the ground 
state. Therefore, we determine T = mdV^/2 + mivf /2 as the average kinetic energy of 
diquark and light quark, and T/2 = m c i^i/2 + m&v^/ '2 as the average kinetic energy of 
heavy quarks inside the diquark (the coefficient 1/2 takes into account the antisymmetry 
of color wave function for the diquark). Finally, we have the following expression for the 
matrix element of the heavy quark kinetic energy: 



We use the value T ~ 0.4 GeV, which results in v\ = 0.146, where the dominant 
contribution comes from the motion of heavy quarks inside the diquark. 

Now we would like to estimate the matrix element of chromomagnetic operator, cor- 
responding to the interaction of heavy quarks with the chromomagnetic field of the light 
quark. For this purpose, we will use the following definitions: O mag = Ya=i 3^r<^ l ° '^G^ 11 
and O mag ~ A(j(j + 1) — Sdi^d + 1) — sj( s Z + where Sd is the diquark spin (as was 
noticed by the authors earlier Q, there is only the vector state of the cc-diquark in the 
ground state of such baryons), si is the light quark spin and j is the total spin of the 
baryon. Since both c-quarks additively contribute to the total decay width of baryons, 
we can use the diquark picture and substitute for the sum of c-quark spins the diquark 
spin. This leads to the parameterization for O mag , as it is given above, and, moreover, it 
allows us to relate the value of the matrix element for this operator to the mass difference 
between the excited and ground state of baryons: 



The account for the interaction of heavy quarks inside the diquark leads to the following 
expressions for the chromomagnetic and "Darwin" terms: 




(21) 




(22) 



(5& ) |fr+ g q-£*& c |5& ) ) 
2M =W 



9 m c 



2 o|*(0)| 2 



(23) 
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1 - 5 2 |*(0)| 2 + ... 



« 1 - 0.074 - 0.004 - 0.003 - 0.005 + . . . (25) 

We can see that the largest contribution to the decrease of the decay width comes 
from the time dilation, connected to the motion of heavy quarks inside the baryon. For 
the matrix element of the operator cga^G^c, we get: 



(^W^clEff) = 4 ^)-^)) 
2M„(*) 3 m c 



'=(*) 

'— 'CC 

An 2 a 2 

Now let us continue with the calculation of the matrix elements for the four-quark op- 
erators, corresponding to the effects of Pauli interference and weak scattering. The 
straightforward calculation in the framework of nonrelativistic QCD gives 

(c7m(1 " 75)c)(g 7 M (l - 75)<?) = 2m c y- 1 (l - 4S c S q ), (27) 

(c7^75c)(97 M (l - 75)9) = -*S c S q ■ 2m c V-\ (28) 

where V~ l = | v I / i(0)| 2 , and ^i(0) is the light quark wave function at the origin of two 
c-quarks. We suppose, that |^i(0)| has the same value as that in the D-meson. So, we 
find: 

|«,<0)|' - (29) 

Then, again remembering that both c-quarks additively contribute to the total decay 
width and using the diquark picture, we can substitute for S C1 + 5 C2 the Sd, where Sd is 
the diquark spin. Thus, we have 

<3W 1(5^(1- 75)c)(*7"(l -7b)9)|sW) = l m c • |^(0)| 2 , (30) 
{^\(c W c)(q^(l- 75 )q)\E c ^) = 8m c • |^(0)| 2 . (31) 

The color antisymmetry of the baryon wave function results in relations between the 
matrix elements of operators with the different sums over the color indexes: 

(^\(- Cl T,c k )(q k r(l ~ 75)%|4* ) > = -(SW|(CV)(^(1 " 75)<7|£W), 

where T M is any spinor structure. Thus, we completely derive the expressions for the 
evaluation of the required matrix elements. 



4 Numerical estimates 

Now we are ready to collect the contributions, described above, and to estimate the total 
lifetimes of baryons E++ and £+. For the beginning, we list the values of parameters, 
which we have used in our calculations, and give some comments on their choice. 

m c = 1.6 GeV, m s = 0.45 GeV, \V CS \ = 0.9745, 

M w++ = 3.56 GeV, M =+ = 3.56 GeV, M* M - M„ w = 0.1 GeV, 

"cc "cc ~} c< J = cc 

T = 0.4 GeV, |^(0)| = 0.17 GeV*, mi = 0.30 GeV. 
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For the parameters M~++, M~+ and M* M — M„u\ we use the mean values, given in 

— cc — cc =cc 

the literature. Their evaluation has been also performed by the authors in the potential 
model for the doubly charmed baryons with the Buchmiiller-Tye potential, and also in 
papers [17, [l^, [l^, |2(J. For fjj we use the value, given in [||, |2lJ and for T we take it from 



[22]. The mass m c corresponds to the pole mass of the c-quark. For its determination we 
have used a fit of theoretical predictions for the lifetimes and semileptonic width of the 
D°-meson from the experimental data. This choice of c-quark mass seems effectively to 
include unknown contributions of higher orders in perturbative QCD to the total decay 
width of baryons under consideration. 

The renormalization scale \i is chosen in the following way: [x\ = m c in the estimate of 
Wilson coefficients C for the effective four-fermion weak lagrangian with the c-quarks at 
low energies and \ii = 1.2 GeV for the Pauli interference and weak scattering (/c-factor). 
The latter value of renormalization scale has been obtained from the fit of theoretical 
predictions for the lifetimes differences of baryons A c , H+, S° over the experimental 
data. Here we would like to note, that the theoretical approximations used in M include 
the effect of logarithmic renormalization and do not take into account the mass effects, 
related to the s-quark in the final state. For the corresponding contributions to the decay 
widths of baryons with the different quark content we have: 

Ar ni (A c ) = c d (O d ) Ac + c u (O u ) Ac , 

Ar nl (~t) = c s (O s ) Et +c u (O u ) E +, (32) 
AT nl (E° c ) = c d (O d ) s o+c s (O s ) 3 o, 

where (O q )x c = {X c \O q \X c ), O q = (07^0) ((fy^t?) and q = u,d,s. The coefficients c g (/i) 
are equal to: 



c d 



G)ml 



[Cl + Cl + i(4fc£-l)(C?-C£)], 
[(C+ + C__) 2 + ~(1 - 4^)(5C 2 + Cl - 6C + C_)], (33) 



4vr L + "3 



16vr '3 

'" lc [(C+ - C_) 2 + \(1 - Ak^)(5Cl + C 2 _ + 6C+C_)] 



G 2 f m 2 



16tt LV ^ 3 

We use the spin averaged value of the D-meson mass for the estimation of the effective 
light quark mass mi as it stands below: 

mD = mc + m + „ 1.98 GeV. (34) 

m c + mi 

The s-quark mass can be written down as: 

m s = mi + 0.15 GeV. (35) 
As we have already mentioned, the spectator decay width of c-quark T C)Spec is known in 



the next-to-leading order of perturbative QCD [12, 13, 14, 15, 16 1. The most complete 



calculation, including the mass effects, connected to the s-quark in the final state, is 



given in [16]. In the present work we have used the latter result for the calculation of 
the spectator contribution to the total decay width of doubly charmed baryons. In the 
calculation of the semileptonic decay width, we neglect the electron and muon masses in 
the final state. Moreover, we neglect the T-lepton mode. 
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Now, let us proceed with the numerical analysis of contributions by the different 
decay modes into the total decay width. In table 1 we have listed the results for the 
fixed values of parameters, described above. From this table one can see the significance 
of effects caused by both the Pauli interference and the weak scattering in the decays 
of doubly charmed baryons. The Pauli interference gives the negative correction about 
36% for the H+ + -baryons, and the weak scattering increases the total width by 144% for 
£+. As it has been already noted in the Introduction, these effects take place differently 
in the baryons, and, thus, they enhance the difference of lifetimes for these hadrons. 

It is worth here to recall that the lifetime difference of D + and L>°-mesons is generally 
explained by the Pauli interference of c-quark decay products with the antiquark in the 
initial state, while in the current consideration, we see the dominant contribution of weak 
scattering. This could not be surprising, because under a more detailed consideration we 
will find, that the formula for the Pauli interference operator for the L>-meson coincides 
with that for the weak scattering in the case of baryons, containing, at least, a single 
c-quark. 

Finally, collecting the different contributions for the total lifetimes of doubly charmed 
baryons, we obtain the following values: 

t = ++ =0.43ps, r = + =0.12ps. 

'CC — 'cc 

Rather broad variations of both the c-quark mass in the range of 1.6 — 1.65 GeV and 
the mass difference for the strange and ordinary light quarks in (35) in the range of 
0.15 — 0.2 GeV, lead to the uncertainties in the lifetimes: 5t~++ = ±0.1 ps, 5t~+ = 

' — 'CC 1 — 'CC 

±0.01 ps. 

5 Conclusion 

In this work we have performed a detailed investigation on the lifetimes of doubly 
charmed baryons E+ on the basis of the operator product expansion for the transi- 
tion currents. For the first time, we have presented the formulae with the simultaneous 
account for both the mass effects and low-energy logarithmic renormalization for the con- 
tributions to the total decay width of baryons, containing heavy quarks, as it is caused 
by the effects of Pauli interference and weak scattering. The usage of the diquark picture 
has allowed us to evaluate the matrix elements of operators derived. Further, we have 
discussed the procedure to choose the values of parameters for the total lifetimes of these 
baryons. The obtained results show the significant role of both the Pauli interference 
and the weak scattering. 

In conclusion, the authors would like to express their gratitude to Prof. S.S.Gershtein 
for usefull discussions and, especially, to Prof. M.B.Voloshin for clear explanations to 
some questions appeared during a walking along this work. We want also to thank for 
the hospitality of DESY-Theory group, where a part of this work was done. 

This work is in part supported by the Russian Foundation for Basic Research, grants 
96-02-18216 and 96-15-96575. 
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Appendix 

In this appendix we present the explicit formulae for the spectator decay of the c- 
quark in next-to-leading order of the perturbation theory with the account for the mass 
effects, related to the s-quark in the final state. 

The coefficients C+ and C_ in the effective lagrangian with the account for the next- 
to-leading order in perturbative QCD acquire the additional multiplicative factors: 

7tt ( 1 1\ 1 , a s{m w ) - ajjj) 7± 0) / 7± } 01 \ , a s (m w ) 

F±(/i) = 1 + T, Wo 
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where 7^ is the coefficients of anomalous dimensions for the operators 0±: 

(0)(x s (i) { a s \ 2 . 3 . 



with 



in the naive dimensional regularization (NDR) with the anticommutating 75, and 
is a number of flavours taken into account, is the initial two coefficients of QCD 
/^-function, 

P = -9s{Po^ + Pi (£) 2 + 0(«2)}, 

2 38 
Aj = 11 - -n/, A = 102 - — n/ . 

The coefficients i?± are written down in accordance to the requirement of agreement 
between the effective lagrangian, evaluated at the scale fi = mw, and the Standard 
Model one up to terms of the order of a 2 (mw)'- 

2jV c 

where N c = 3 is the number of colors. In the NDR scheme for B, we find B = 11. 

Using the effective lagrangian in the next-to-leading order of perturbative QCD and 
calculating the one-gluon corrections, we get the following expression for the spectator 
c-quark decay: 

T(c^sud) = r [2C^(/x) + C 2 (fi) + a ^ mw l ~ a ^ {2C 2 + ( f i)R + + + 



2vr 



-{2C^(/i)5 + + Cl(/i) J B_} + 



where 



- C^)} 22 -^{G C + G d } + (36) 
\{C 2 + {„) - C 2 + ^)} 2 -^{G a + G b + Ge}}, 



3 „4 10 2 



/i(o) = l-8a + 8a J -a 4 -12a^lna, 



and 



7f /t£ A 
2/3o I 7 £°) A) J ' 
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For G a ,Gb,G c , Gd and G e we have found: 

{G a + G h )h{a) = ^ -7T 2 -a[80-lno] + 32a 3/ V 

-a 2 [273 + 16tt 2 - 18 In a + 36 In 2 a] + 32a 5/ V 
-^a 3 [118-571na]+0(a 7/2 ), (37) 

{G c + G d )h{a) = ^-7r 2 -8a[10-7r 2 + 31na]-a 2 [117 - 247r 2 + 

(30 - 8vr) In a + 36 In 2 a] - ^a 3 [79 + 2tt 2 
-62 In a + 6 In 2 a] + 0(a 4 ), (38) 
{G a + G b + G e + B)h{a) = (6ln^ + ll)f 1 (a)-^--7T 2 + 8a[21-TT 2 -3lna] 

+32a 3/ V - a 2 [111 + 40vr 2 - 258 In a + 36 In 2 a] 

+32a 5/ V - ^a 3 [305 + 18vr 2 + 30 In a - 54 In 2 a] 
9 

+0(a 7/2 ). (39) 

These approximations can be used in the range of a values: a < 0.15, where a = (^-) 2 , 
which, indeed, takes place in the calculations under consideration. 
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Mode or decay 


Width, ps" 1 


Contribution in % 


Contribution in % 


mechanism 






(S+) 


Cspec ^ sdll 


2.894 


124 


32 


c — ► se + f 


0.380 


16 


4 


PI 


-1.317 


-56 




WS 


5.254 




59 




2.337 


100 






8.909 




100 



Table 1: The contributions of different modes to the total decay width of doubly charmed 
baryons. 
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Figure 1: The spectator contribution to the total decay width of doubly charmed baryons. 
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Figure 2: The Pauli interference of c-quark decay products with the valence quark in the 
initial state for the S+ + -baryon. 
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Figure 3: The weak scattering of the valence quarks in the initial state for the S+-baryon. 
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